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Abstract. We estimate the diffuse supernova neutrino background (DSNB) using the recent 
progenitor-dependent, long-term supernova simulations from the Basel group and including 
neutrino oscillations at several post-bounce times. Assuming multi- angle matter suppression 
of collective effects during the accretion phase, we find that oscillation effects are dominated 
by the matter-driven MSW resonances, while neutrino- neutrino collective effects contribute at 
the 5-10% level. The impact of the neutrino mass hierarchy, of the time-dependent neutrino 
spectra and of the diverse progenitor star population is 10% or less, small compared to the 
uncertainty of at least 25% of the normalization of the supernova rate. Therefore, assuming 
that the sign of the neutrino mass hierarchy will be determined within the next decade, the 
future detection of the DSNB will deliver approximate information on the MSW-oscillated 
neutrino spectra. With a reliable model for neutrino emission, its detection will be a powerful 
instrument to provide complementary information on the star formation rate and for learning 
about stellar physics. 
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1 Introduction 

The Diffuse Supernova Neutrino Background (DSNB) is the flux of neutrinos and antineutri- 
nos emitted by all core-collapse supernovae (SN) in our universe. Isotropic and stationary, it 
is a guaranteed signal that will give us a unique image of the SN population of the universe 
and of its history all the way to redshift ~ 1 or so. This image can be a precious alternative, 
or a unique complement, to a high statistics detection of an individual nearby SN which is 
a rare event (~ 1-3 SNe per century in our galaxy [1-4]). Although the DSNB has not been 
detected yet, its discovery prospects are excellent. Gadolinium enriched Super-Kamiokande 
is expected to detect the DSNB at a level of few events per year [5, 6] and tens to hundreds of 
events are predicted for the next generation of larger scale, 0.1-1 Mt mass detectors, running 
for a few years [7, 8]. 

The DSNB depends on the cosmological SN rate, the neutrino fluxes and the oscillation 
physics. Its prediction is a challenging problem because of the still sparse astrophysical 
observations of supernovae, the incomplete picture of neutrino oscillation physics, and several 
obstacles — of theoretical and computational nature — in the way of fully modeling the 
neutrino emission from core collapse. Due to these difficulties, a number of simplifying 
assumptions were adopted to model the DSNB, typically the same neutrino spectra were 
used for different post-bounce times and for all progenitors and only approximate neutrino 
oscillation physics was considered (see [9, 10] for reviews on the topic). Only recently efforts 
were developed to include more accurate astrophysical information and oscillations driven by 
neutrino- neutrino forward scattering [11-14]. 

From the perspective of neutrino oscillation physics, lately new advances have taken 
place on the experimental and theoretical front. The elusive mixing angle #13 has been 
measured [15, 16], leaving us only with open questions on the neutrino mass hierarchy and 
the CP-violating phase (although the latter does not affect neutrino oscillations in super- 
novae [17]). Concerning neutrino oscillations in supernovae from the theoretical perspec- 
tive, the matter enhanced flavor conversion picture [18, 19] has been completed by includ- 
ing neutrino-neutrino interactions, responsible for the non- linear neutrino flavor evolution 
(see [20] for a review on the topic). The present sketch is that, during the accretion phase, 
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total (or partial) multi-angle matter suppression of collective oscillations is expected [21-26]; 
while, during the cooling phase, multiple spectral splits occur [27, 28] and Mikheev-Smirnov- 
Wolfenstein (MSW) resonances are expected at larger radii [18, 19]. Moreover, most recent 
numerical simulations of core-collapse SN provide information about neutrino spectra up to 
~ 10 s post-bounce [29, 30] and for different progenitor masses [30]. Therefore the time is 
mature to include these updates in the estimate of the DSNB and to assess the impact of 
each in comparison with the many uncertainties, experimental and theoretical, that affect the 
diffuse flux itself. This is the goal of this paper. We model the DSNB, as consistently as pos- 
sible, in the framework of the Basel model [30]. We consider oscillation physics for different 
post-bounce times for both neutrino mass hierarchies, including the multi-angle treatment of 
v-v interactions. 

This work is organized as follows. In Section 2 we discuss the cosmological super- 
nova rate, the adopted supernova models for different progenitor masses and we give a brief 
overview on neutrino oscillation physics in supernovae. In Section 3 we discuss the impact 
of neutrino oscillations on the neutrino time- integrated flux for a fixed supernova mass. In 
Section 4, we present our results on the DSNB. Conclusions and perspectives are illustrated 
in Section 5. 

2 Input physics 

In this Section, the cosmological SN rate and the reference neutrino signal for different 
SN masses are discussed. We also introduce the neutrino mass-mixing parameters and the 
quantum kinetic equations describing the neutrino oscillations in supernovae. 

2.1 Diffuse SN neutrino flux 

Simply put, the DSNB offers us an image of the entire SN population of the universe. It 
reflects the demographics of this population: how it is distributed in space, its history, the 
diverse sub-types that contribute to it with their different neutrino emissions. Besides, the 
image we get is influenced by the neutrinos' own history, as they travel cosmic distances 
before they reach us. Therefore the diffuse flux of a given neutrino flavor depends on many 
physical quantities of different origin. 

• The flux produced by an individual SN, which depends on the dynamics of the core 
collapse and is different for different masses of the progenitor star. 

• Propagation effects: namely flavor oscillations, which depend on the neutrino masses 
and mixings, and the redshift of energy; for a neutrino produced at redshift z with 
energy E', the observed energy is E = E' / (1 + z). 

• The comoving SN rate (SNR), psn^, M), per unit of redshift and per unit of progen- 
itor mass M. It is defined in the mass interval between Mq ~ 8 M Q (the minimum 
necessary to have a core-collapse supernova) and a maximum M max ~ 125 M@, the ten- 
tative upper limit for the occurrence of normal core-collapse supernovae (as opposed 
to pair instability ones or black- hole forming events) [11]. In terms of redshift, the 
SNR is largest between z = and z max — 5, reflecting the period of most intense star 
formation activity (see Sec. 2.2). Note that we neglect the contribution due to failed 
supernovae [31]. 
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With these ingredients, the diffuse flux for each flavor ug (/3 = e, e, p, or r) can be written 
as [11]: 

_ c f M ^ /■*-« Psn(z, M)F Vj} (E', M) 

n oJm Jo V"m(1 + 2) d + "A 

where c is the speed of light and i^o = 70 km s _1 Mpc -1 the Hubble constant; Qm = 0.3 
and Oa = 0.7 are the fractions of the cosmic energy density in matter and dark energy 
respectively and F U „(E' , M) is the oscillated ug flux for a SN with progenitor mass M (see 
Sec. 3). As we will see, the DSNB is dominated by the contribution of the closest (z 1) 
and least massive (M ~ Mo) stars, and depends only weakly on M max and 2 max . 

Although the DSNB has not been observed yet, interesting upper limits exist. The 
most stringent is on the v e component of the flux, from a search of inverse beta-decay 
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events at Super-Kamiokande, above 17.3 MeV threshold: (j) Pe & 2.8-3.0 cm s at 
C.L. [32, 33]. This bound is generally consistent with predictions, excluding scenarios where 
multiple parameters conspire to generate a particularly large flux [34]. Because the search at 
Super-Kamiokande is background-dominated, any substantial improvement on it will require 
better background subtraction. Methods involving water with Gadolinium addition [5], liq- 
uid Argon [35-37], and liquid scintillator [38] are especially promising. Of these, detectors of 
Megaton class will have the further advantage of high statistics, yielding up to hundreds of 
events a year from the DSNB (see e.g., [7, 8]). 

2.2 Cosmological supernova rate 

Considering that SN progenitors are very short lived, the SNR is proportional to the Star 
Formation Rate (SFR), p+, defined as the mass that forms stars per unit time per unit 
volume. The relationship between the SNR and SFR is given by the Initial Mass Function 
(IMF), <q(M) oc M -2 - 35 [39], which describes the mass distribution of stars at birth: 

PSN (z, M) = Uz) . (2.2) 

CmI dM M V (M) 

Recent analyses of SNR and SFR measurements [40] show that a piecewise parametrization 
of p\ is adequate [41]: 

({l + z) s z<l 
p± oc I [l + z) a 1< z < 4.5 . (2.3) 
I (1 + z)7 4.5 < z 

Here we adopt this function, with S, a, 7 and the normalization fixed at the best fit values [41] : 
S = 3.28, a = -0.26, 7 = -7.8 and J^ max dM psn(0, M) = 1.5x 10~ 4 Mpc^yr- 1 (integrated 
SNR at the present epoch). Note that, due to the redshift of energy, neutrinos of higher 
redshift accumulate at lower energies, so that in the energy window relevant for experiments 
(11 MeV ^ E ^ 40 MeV [6, 32]) the diffuse flux is dominated by the low z contribution, 
z ^5 1. Therefore, its dependence on a and 7 is weak. The flux is also dominated by the 
lower mass stars, considering the fast decline of the IMF with M, so that there is a strong 
dependence on Mq, but a weak one on the high cutoff M max . 

Let us now comment on the existing measurements of the SNR and their uncertainties, 
which are dominated by normalization errors [3, 40]. Perhaps the most precise way to measure 
the SNR is from data on the SFR, via Eq. (2.2). The cosmic star formation history as a 
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function of the redshift is known from data in the ultraviolet and far-infrared [41]. Constraints 
are stronger at z ^ 1, while observations at higher redshift have larger errors. Following 
Ref. [40], we adopt the statistical error of about 25% as uncertainty on the normalization. 
This underestimates the total error, that includes a number of systematic effects, some of 
which are difficult to estimate. 

In principle, to rely on the SFR data is not necessary: the SNR is given by direct SN 
observations. However, data from SN surveys are sparse (low statistics) and extend only to 
z ~ 1. The most precise measurement is the recent LOSS result for the local (z = 0) rate, 
integrated over all masses [42]: J^ max dM p SN (0,M) = (0.70 ± 0.15) x 10~ 4 Mpc~ 3 yr- 1 . 
Here we use its error (dominated by systematics), of ~ 22% as a very optimistic estimate 
of uncertainty. Other errors are due to possible dust obscuration and supernova impostors 
(see [40] for details). 

Surprisingly, the normalization from direct SN observations is lower than that from SFR 
data, by a factor ~ 2 (with a ~ 2a level significance at z = 0) and by a smaller factor at 
higher z [3]. A possible explanation to this mismatch is that many supernovae are missed 
because they are either optically dim (low-luminosity) or dark, whether intrinsically or due 
to obscuration. Other proposed explanations are an incomplete understanding of the star 
formation and supernova rates including that supernovae form differently in small galaxies 
than in normal galaxies [3]. 

Finally, a complementary constraint on the SNR normalization comes from the Super- 
Kamiokande limit on the DSNB (Sec. 2.1). The constraint depends on the neutrino spectrum 
and is strong for rather hot neutrino spectra [43]. For the cooler spectra favored by SN 1987A, 
data allow a normalization up to ~ 2 times higher that the value we use here. 

2.3 Reference neutrino signals from supernovae of different masses 

We adopt supernova simulations consistently developed over 7 s or more for three different 
masses (8.8 M , 10.8 M and 18 M & ) [30]. Such core-collapse models are based on spher- 
ically symmetric general relativistic hydrodynamics and they include spectral three-flavor 
Boltzmann neutrino transport with the equation of state from Shen et al. [44]. However, 
they do not include nucleon recoil. As observed in [45], the inclusion of nucleon recoil is 
responsible for reducing the differences among the mean energies of different flavors during 
the cooling phase (otherwise different up to 10%) [29, 45, 46]. Since one of our purposes 
is to evaluate the impact of neutrino oscillations on the DSNB, the Basel model favors the 
appearance of any signature due to neutrino oscillations in the DNSB . 

At a radius r, the unoscillated spectral number fluxes for flavor up and for each post- 
bounce time i p b are 

o _ L up (t ph ) tp vp (E,tp h ) _ F^(E,t ph ) 

where L„_(ipt,) is the luminosity for flavor up, {E u ^(t p b)) the mean energy, and Lp V/3 {E,t p \ > ) a 
quasi-thermal spectrum. We describe it schematically in the form [47] 

^(E,t ph ) = fcfc b ) {jE^) MtVh) e-MWW&'W • (2.5) 

The parameter ap(t ph ) is defined by (E Ufi (t ph ) 2 } / ' {E Up (t ph )) 2 = (2 + ap(t ph ))/(l + ap(t ph )) 
and Cp(t p b) is a normalization factor such that f dE tp U/3 (E,t p b) = 1- 
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Figure 1. Luminosities (on the top) and mean energies (on the bottom) in the observer frame for 
three progenitors with masses M = 18, 10.8, 8.8 Mq (from left to right) as a function of the post- 
bounce time [30, 48]. In blue (red, black respectively) are plotted the quantities related to v^ T and 
v^.r denoted by v x (v e , v e ). 

Figure 1 shows the luminosities (on the top) and the mean energies (on the bottom) in 
the observer frame for the three adopted progenitor models as a function of the post-bounce 
time as in [30, 48]. Note that the luminosities of the different flavors are almost equal during 
the accretion phase, while L Ve 3> Ly^L v during the cooling phase. While, the mean 
energies are (E v ) > (Ep e ) > (E Ue ) during the cooling phase. 

2.4 Neutrino mixing parameters and quantum kinetic equations 

We assume the following neutrino mass squared differences [49] 

6ml tm = 2.35 x 1(T 3 eV 2 , (2.6) 
<5m 2 ol = 7.58 x 1(T 5 eV 2 , (2.7) 

and we discuss both normal (NH, <5m 2 tm > 0) and inverted hierarchy (IH, £m 2 tm < 0) 
scenarios. The mixing angles are [49, 50] 

sin 2 6»i3 = 0.02 and sin 2 12 = 0.3 ; (2.8) 

we neglect the third mixing angle 623 for reasons that will be clear in a while. 

We treat neutrino oscillations in terms of the matrices of neutrino densities for each 
neutrino mode with energy E, pg, where diagonal elements are neutrino densities, off-diagonal 
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elements encode phase information due to flavor oscillations. The radial flavor variation of 
the quasi-stationary neutrino flux is given by the "Schrodinger equation" 

id r pE = [He,Pe] and id r p E = [H#, p E ] , (2.9) 

where an overbar refers to antineutrinos and sans-serif letters denote 3x3 matrices in flavor 
space consisting of u e , and v T . The initial conditions are p E = diag(n I/e , n u , n Vr ) and 
Pe = diag(n Pe , r2,p M ,np r ) where n Uli — tl Ut — n Ux and the same for antineutrinos. The 
Hamiltonian matrix contains vacuum, matter, and neutrino-neutrino terms 

\\ e = WT + Ue + Ue ■ (2-10) 

In the flavor basis, the vacuum term is a function of the mixing angles and the mass-squared 
differences 

Hf c = Udiag(-^,+^ H )lJt , (2.11) 

where U is the unitary mixing matrix, function of the mixing angles, transforming between 
the mass and the interaction basis, ujl = 5m^ oi /2E and ujh = ^ m ttra/^- The matter term 
includes charged-current (CC) interactions and it is responsible for MSW resonances. In the 
flavor basis, it is 

H m = V2G F diag(iV e , 0, 0) , (2.12) 

where N e is the net electron number density (electrons minus positrons). 

The corresponding 3x3 matrix caused by neutrino- neutrino interactions [51] including 
the multi-angle dependence is 

= f dE'( PE i - PE>){1 - COS0EE') , (2.13) 

where Qee 1 is the angle between the momenta of the colliding neutrinos [52]. 

Ignoring radiative corrections, Un and v T are exactly equivalent in supernovae, allowing 
us to define an interaction basis (u e , v x , u y ) = R T (92z) {v e , u T ), with i? T (^23) the transpose 
of the rotation matrix 11(623), such that effectively #23 = [53]. Note that, in this rotated 
basis, the e — y oscillations are driven by 6m^ tm while the e — x oscillations are driven by 

^sol- 
It is possible to factorize the effects of self-induced flavor conversions and the MSW 

resonances in most of the cases since they are occurring in well separated regions. We 

consider multi-angle v-v interactions driven by the atmospheric mass difference <5m a t m and 

by the mixing angle #13 between v e and u y , while the other flavor (y x ) does not evolve. In 

fact the effects on collective flavor conversions induced by the third flavor are mainly of the 

nature of a subtle correction [20, 54-56] negligible for our purposes. The only way v x could 

affect the final neutrino spectra is by MSW transitions, occurring at larger radii than the 

ones where collective effects happen. Thus the fluxes after the collective oscillations (F^) 

are given by 

= P c Fl + (1 - P C )F° and F£ = P C F° + (1 - P c )F° y , (2.14) 

where P c and P c are the v e and v e survival probabilities after the multi-angle self-induced 
flavor conversions. They are functions of the neutrino energy, and depend on the mass hier- 
archy. The oscillated v y and v y fluxes, F£ and F£ y , can be estimated from the conservation 
of the total number of v (and v): 

Fl y + Fl = F° + F° and i% + i% = F° + F° y . (2.15) 
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The third state, v x (v x ), is not affected by the collective evolution, therefore ,_ ^ = v 
As we consider the self-induced neutrino oscillations as factorized from the MSW in 
first approximation, the fluxes F£ will undergo the traditional MSW conversions after v—v 
interactions. In NH, the MSW resonance due to <5m 2 tm affects the v e flux while the v e flux 
remains almost unaffected. On the other hand for IH, the same resonance affects the v e flux 
and not the v e flux. The fluxes (F Ue and Fp e ) reaching the earth after both the collective and 
MSW oscillations for NH and IH and for large #13 are [13, 53, 57]: 



F™ = sm 2 e l2 [l-Pc(F^Fi,E)}(F» e -F Uy ) + F° Uy , (2.16) 

F™ = cos 2 e 12 P c (F^,Fl,E)(Fl - F° y ) + F° , (2.17) 

F™ = sin 2 9 12 P c (FS e , Fl, E){Fl - F°J + F° , (2.18) 

= cos 2 ^ 12 [l- J P c (^ e ,^,F)](F° + . (2.19) 



Here we have used the fact that, by combining Eqs. (2.14, 2.15) with J dE(F® — Fp e ) = 
const., one can express P c and P c , for each energy, as functions of the fluxes after collective 
oscillations: P c (F£ e , F£ e , E) and P c (F£ e , F£ e , E). These probabilities exhibit a well known 
step-like behavior, that appears in the fluxes F Ve and Fp e as the so called "spectral splits" [28, 
58]. In reality, a number of effects smooth out the splits in the observed neutrino fluxes; we 
discuss this point further below. 

3 Time-integrated neutrino fluxes and oscillation effects 

As a first step towards the calculation of the DSNB, we compute the oscillated fluxes, Eqs. 
(2.16)-(2.19), for each progenitor model, at fixed time snapshots (as an approximation of the 
continuous time evolution, which is too demanding for state of the art computers). The results 
are then used to obtain the up and v$ fluxes integrated over the duration of the neutrino 
burst. For illustration, in this section we discuss the results for the 10.8 Mq progenitor; in 
Sec. 4 the fluxes for all progenitors will be summed up to obtain the DSNB, via Eq. (2.1). 

While the MSW effect is well described analytically, the collective effects require a 
numerical calculation to obtain the probabilities P c and P c . Let us discuss them here in more 
detail. 

The spectral split patterns (affecting P c and P c ) are known to be crucially dependent on 
the initial relative flux densities and on the mass hierarchy. For definiteness, it is convenient 
to distinguish between the probabilities in the accretion phase (i p b < 1 s), and those in the 
cooling phase (i p b ^ Is). In the accretion phase, the multi-angle effects associated with dense 
ordinary matter suppress collective effects [23-26]. Therefore, we adopt the results presented 
in [23, 24] for the two-flavor system (u e ,u y ), considering partial or no flavor conversion for 
several t p y, as in [23, 24]. 

During the cooling phase, the fluxes of different flavors are slightly different, and spectral 
splits occur for neutrinos and/or antineutrinos according to the mass hierarchy and to the 
number of crossings in the non-oscillated spectra (i.e., energies where F£(F) = F® x (E) and 
the same for v) [27]. We calculate these effects by numerically solving Eqs. (2.9) for t p ^ = 
1, 3, 6, 9 s, for the system {y e , v y ). Concerning the neutrino emission geometry, we assume a 
spherically symmetric source emitting neutrinos and antineutrinos like a blackbody surface, 
from a neutrinosphere with radius, R u , that varies with t p b- We define the neutrinosphere 
radius as the radius where the neutrino radiation field is half-isotropic [24, 60] and we adopt 
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Figure 2. Time-integrated E 2 x F Vp (E) fluxes as a function of the energy for M = 10.8 M®, for v e 
(v e ) on the left (right) and NH (IH) on the top(bottom). The non-oscillated fluxes are plotted with 
dotted lines, the fluxes obtained including neutrino interactions with matter only (MSW) are plotted 
with dashed lines and the fluxes obtained including both matter and v background (MSW+^-^) are 
plotted with solid lines. 



the same R u for different flavors. Specifically, we use Sm^ tm as in Eq. (2.6), R u ~ 13 km 
(although slightly varying with i p b) and, to be coherent with the treatment in [23, 24], we 
simulate the matter effect on collective oscillations using a small, effective, in-medium mixing 
angle 6is e S = 10 -3 [61]. The non-oscillated spectra are defined as in Eq. (2.5) with mean 
energies and luminosities as in Fig. 1. 

We find that that multiple crossings appear in the non-oscillated spectra and there- 
fore multiple spectral splits are expected during the cooling phase due to collective effects. 
However, the shape of the splits is smeared and their size is reduced due to the similarity 
of the unoscillated flavor spectra and to multi-angle effects. Often, only partial conversion 
is realized. In particular, for NH we find that P c exhibits a transition from high (P c 0.6) 
to low {P c £ 0.3) values as the energy increases, with transition energy increasing with time 
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from ~ 15 to ~ 35 MeV. A similar behavior is observed for P c for energies of relevance for 
detection, E^,5 — 7 MeV. For IH, P c has a somewhat opposite energy dependence, changing 
from P c ~ 0.1 — 03 at E ^ 15 — 20 MeV to P c ^ 0.6 at higher energy. The transition becomes 
gentler in energy as time increases. For antineutrinos, we have P c ^ 0.7 at all times and all 
energies, with only a moderate energy dependence. Our results for P c , P c are qualitatively 
similar to what shown in previous literature, although smoothened by the multi-angle depen- 
dence and depending on the adopted supernova progenitor; therefore we choose to be brief 
here, and refer to dedicated papers for details [28, 62]. 

To obtain the time-integrated neutrino fluxes, we approximate the time dependence of 
P c and P c with a step-like form, taking them to be constant in time intervals centered on the 
scanned i p b (we adopt a time window of ~ 3 s centered on each i p b). We checked that the 
particular choice of the time-binning affects the final results only weakly. The time-integrated 
fluxes have been obtained by integrating over time the fluxes computed at each time slice. 

Figure 2 shows the time-integrated fluxes, E 2 x F U/3 (E), as a function of the energy for 
v e (z/ e ) on the left (right), for both mass hierarchies for the 10.8 Mq supernova model. We 
distinguish among non-oscillated integrated fluxes (dotted lines), oscillated fluxes including 
MSW effects only (dashed lines) and oscillated fluxes including both MSW and v-v interac- 
tions (solid lines). Note that the MSW effect is large and generates hotter fluxes. Compared 
to MSW only, the collective effects produce a slight hardening or softening of the fluxes de- 
pending on the hierarchy, according to Eqs. (2.16-2.19). No signature of multiple spectral 
splits appear in the time-integrated flux, because they are suppressed during the accretion 
and they occur at different energies for each i p b during the cooling; moreover multi-angle 
effects smear the conversion probabilities. 

Overall, the flux variation due to collective effects - relative to MSW conversion only 
- is at the level of 10% or less. Such a small difference is explained by the fact that, in 
contrast with earlier approximate results on collective effects, the probabilities P c and P c 
rarely approach zero, but remain in the interval ~ 0.3 — 0.8 in most cases. One should also 
consider that the contribution of P c and P c appears in Eqs. (2.16-2.19) weighed by factors 
of sin 2 9i2 or cos 2 #12, that reduce its size. A further reduction takes place in the time- 
integrated flux, which receives about 1/2 of its luminosity from the accretion phase, where 
collective effects are suppressed. Finally, the differences between the non-oscillated fluxes in 
the different flavors become progressively smaller in the cooling phase, thus suppressing any 
spectral change due to oscillations. 

To further illustrate the size of each effect we included, in Tab. 1 we provide the mean 
energies and the mean squared energies for the time-integrated spectra. In the existing 
literature, one progenitor has been often considered as representative of the whole stellar 
population and one post-bounce time was taken as representative of both accretion and 
cooling phase. Therefore, for sake of completeness, we give this case (case a in the Table) 
using the parameters for i p b = 0.5 s as representative of all post-bounce times and the 
10.8 Mq progenitor as representative of the whole stellar population. Looking at the average 
energies in the Table, we see an expected decrease (5-10% or so) when the full cooling over 
~ 10 s is included (case b), and a marked hardening due to including the MSW conversion 
(up to 30%, case c). As expected from Fig. 2, adding collective effects amounts to a change 
of 10% or less (case d), and only a very minor hardening is due to summing the results 
for different progenitor masses (case e). In the next Section, we will discuss in detail the 
differences among the different cases for the DSNB. 
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162.56 


170.62 


174.48 


(Eo e , 


>IH 


(MeV) 


11.82 


10.83 


12.24 


12.17 


12.18 


(ei: 


>IH 


(MeV 2 ) 


174.20 


149.17 


209.53 


206.77 


206.06 



Table 1. Mean energies and mean squared energies for the the time-integrated F vp (E). The different 
cases considered are: (a) F Vfi (E) obtained considering £ p b = 0.5 s as representative of all post-bounce 
times and the 10.8 Mq progenitor as representative of the whole stellar population, no oscillations; 

(b) F v „ (E) obtained considering time-dependent fluxes for the 10.8 Mq model, no flavor oscillations; 

(c) F V p(E) obtained considering time-dependent fluxes for the 10.8 Mq model, MSW flavor conver- 
sions; (d) F vp (E) obtained considering time-dependent fluxes for the 10.8 Mq model, MSW + v—v 
interactions; (e) F Vf} (E) obtained including the whole stellar population, time-dependent fluxes, MSW 
+ v-v interactions. 

4 Diffuse supernova neutrino background 

After computing the time-integrated spectra for each of the considered SN progenitors as 
explained in Sec. 3, we derive the DSNB, using Eq. (2.1). We adopt the 8.8 Mq SN as 
representative of supernovae with mass in the interval [8, 10] Mq, the 10.8 Mq SN for masses 
in the interval [10, 15] Mq and the 18 Mq progenitor for M > 15M . Note that the for the 
first interval we adopt the widest range expected for O-Ne-Mg core formation. 

Figure 3 shows the DSNB, E 2 x <& Vg with § Vf} , as a function of the energy for v e (u e ) 
on the top (bottom). In the plots on the right, we focus on the region of interest for Super- 
Kamiokande detection (E > 17.3 MeV) and it is possible to appreciate a maximum variation 
of 10-20% (at E ~ 20 MeV) related to the mass hierarchy. In order to favor a comparison 
with the existing literature, we show the fluxes obtained by taking the 10.8 Mq progenitor as 
representative of all progenitors (thin curves), and the full result that includes the variation 
of the fluxes and probabilities with the progenitor mass (thick lines) . Due to the similarity of 
the initial luminosities and mean energies (see Sec. 3) the differences between the two cases 
are of few per cent. 

Shock effects in Fe-core supernovae manifest themselves during the cooling phase and at 
large radii. We choose to neglect such effects for the 18 Mq and 10.8 Mq models because they 
practically might affect the <5m 2 tm -driven resonance only for the late cooling phase having a 
negligible impact on the DNSB. On the other hand, the 8.8 Mq model has a Ne-O-Mg core 
and the first hundreds of milliseconds of the burst are affected by variations of the matter 
density profile due to the shock passage. The shock wave is responsible for turning the flavor 
conversions from non-adiabatic to adiabatic within the first ~ 250 ms affecting the Sm 2 tm - 
driven resonance [63]. However, a very small time- window would have been affected by the 
shock passage in the electron-capture SN and we expect its impact on the total DSNB to be 
negligible. 

Recently, it has been observed that the non-forward neutrino scattering contribution 
might affect flavor conversions during the accretion phase [64, 65]. These studies reveal that 
the consequences of non-forward scatterings are negligible for large SN masses when com- 
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Figure 3. DSNB as a function of the energy (E 2 x ® Ve ,v^) on the left panels, zoom in the region of 
Super-Kamiokande detection in the right panels (E > 17.3 MeV). The v e (p e ) flux is plotted in black 
(red). The solid (dashed) line represents the DSNB in IH (NH). To mark the effects induced by the 
inclusion of different SN model, the thick lines represent the DSNB including the three different SN 
progenitors and the thin line is the DSNB obtained extending 10.8 M® to the whole stellar population. 

plete matter suppression occurs [65] and for low mass supernovae as the Ne-O-Mg ones [64] . 
However, it might be still relevant for supernovae with intermediate masses (as the 10.8 M 
model). Since these studies are still embryonal and no numerical implementation of the 
non- forward contribution is available yet, a priori, one could expect larger flavor conversions 
induced by the non-forward term during the accretion phase especially for the intermediate 
SN mass range. Therefore, we chose to adopt a parametric study of the accretion phase 
for the 10.8 Mq model, considering two extreme cases of full flavor conversions (P c = 0) 
and complete matter suppression (P c = 1) due to multi-angle matter effects, other than the 
intermediate scenario presented in [23-25]. We find that, at the level of the DSNB, after 
summing the fluxes from progenitor of different masses, the difference between the two cases 
(not shown here) is of the order of few per cent. 

Figure 4 gives results for the DSNB integrated in energy above the experimental thresh- 
old 17.3 MeV [32]. It shows results for the cases (a-e) of Tab. 1 for both mass hierarchies. 
The differences between these cases reflect the results found for an individual SN burst in 
Sec. 3, with more emphasis on the high energy part of the neutrino spectrum (due to our high 
integration threshold) . Compared to static spectra (case a) , the effect of the time-dependence 
of the spectra over the 10 s of the neutrino burst is responsible for a variation of 6% (case b). 
The MSW effects (case c) are the largest source of variation of the DSNB with respect to the 
case without oscillations and it is about 50-60%. This variation is particularly evident for 
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Figure 4. Energy integrated DSNB (for E > 17.3 MeV) as a function of the discussed contributions, 
in black (red) for v e (i> e ). We adopt triangles (dots) for the IH (NH) scenario. From left to right the 
source of variation of the DSNB are: (a) DSNB for t p b = 0.5 s as representative of all post-bounce times 
and the 10.8 Mq progenitor as representative of the whole stellar population, no oscillations; (b) DSNB 
obtained considering time-dependent fluxes for the 10.8 Mq model, no flavor oscillations; (c) DSNB 
obtained considering time-dependent fluxes for the 10.8 Mq model, MSW flavor conversions; (d) 
DSNB obtained considering time-dependent fluxes for the 10.8 Mq model, MSW + v—v interactions; 
(e,f) DSNB obtained including the whole stellar population, time-dependent fluxes, MSW + v—v 
interactions; the difference in the two points is due to the factor of ~ 2 between the results of SN 
surveys and those of star formation measurements (see Sec. 2.2). The error bars are the errors on the 
normalization of the supernova rate [3] . 



neutrinos and the NH, where the difference between the unoscillated and oscillated spectra 
is the largest, due to the complete flavor permutation driven by the large #13. Neutrino- 
neutrino interactions (case d) are responsible for a variation of 5-10% with respect to the 
MSW only case. Summing over the stellar population (case e) is responsible for a DSNB 
variation of 5-10% due to the more luminous fluxes of the more massive stars. For sake of 
completeness, we provide the numerical values for the DSNB for the case (e): 

$ t ^ NH = 0.31 cm-V 1 and $ t ^ IH = 0.27 cnrV 1 , (4.1) 
$ t ^ NH = 0.26 cm-V 1 and $^ IH = 0.32 cm-V 1 . 

The maximum impact given by the mass hierarchy is 20%, triangle and dot in (e). It is 
realized for antineutrinos because, for large #13, the high density MSW resonance is adiabatic 
and the v e flux changes from almost complete survival for NH (survival probability ~ 0.7- 
0.8) to almost complete conversion for IH (survival probability ~ 0.1-0.2). We estimate the 
DSNB taking into account the unexplained factor of ~ 2 between the results of SN surveys 
and those of star formation measurements as discussed in Sec. 2.2, points (e) and (f). This 
mismatch is responsible for the largest astrophysical source of error on the estimation of the 
DSNB, ~ 50%. Moreover, the error on the normalization of the supernova rate is about 25% 
and it is represented by the error bars in (e) and (f). However, this error could be most likely 
higher once several systematic errors are included. 

Another potential source of error on the DSNB is in the equation of state of nuclear 
matter used in core collapse simulation (assumed fixed in all our computations). The variation 
of the total neutrino energy release during the supernova explosion reflects the variation of the 
gravitational binding energy of the neutron star in dependence of different nuclear equations 
of state (see [66] for details). However, we expect that the differences of the time-integrated 
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emitted neutrino spectra using different equations of state will be not large enough to play a 
role in comparison will the other astrophysical uncertainties discussed in this paper. 

5 Conclusions 

The diffuse supernova neutrino background is an integrated picture of all the core collapse su- 
pernovae of different progenitor mass and redshift, and has realistic chances to be detected in 
the near future by upcoming experiments. The predictions for this flux have been constantly 
updated to include new understanding of SN physics and astrophysics, and of neutrino phe- 
nomenology. We present one such update, to include, as consistently as possible, the results 
of new supernova simulations and oscillation calculations. 

Specifically, our calculation is the first that uses neutrino fluxes that are numerically 
calculated over the entire ~ 10 s of the neutrino emission for a set of different progenitor 
masses. For each progenitor star, we study neutrino oscillations at different times post 
bounce, for both the neutrino mass hierarchies, including multi-angle collective effects and 
adopting the recently measured value of #13. The results are then used to obtain the time- 
integrated flux from a single SN; finally, the oscillated fluxes for different progenitor masses 
are combined together into a prediction for the DSNB. 

The results of this paper rely on the Basel simulation of neutrino emission [30]. Al- 
though model dependent, they are indicative of what might be expected from state of the 
art simulations. In particular, the Basel model has the largest spectral differences between 
the neutrino flavors, and therefore it is useful to estimate the largest oscillation effects that 
can be realistically expected. 

Let us summarize our main results. 

• The inclusion of time-dependent neutrino spectra are responsible for a colder neutrino 
spectrum compared to calculations using a fixed spectrum from early time emission 
(e.g., t p b ~ 0.5 s). For a single supernova, the effect is a ~ 5% shift in the average 
energies (Tab. 1), which translates into a ~ 5% difference in the integrated DSNB above 
a 17.3 MeV threshold (Fig. 4). 

• The largest effect of flavor oscillations is due to the matter-driven MSW resonances (~ 
50-60%), while neutrino-neutrino collective effects give a further ~ 5-10% contribution. 
Moreover, the DSNB does not present any energy-dependent signature of the collective 
oscillations. For fixed SN progenitor, multi-angle matter effects suppress collective 
oscillations in the accretion phase, while multiple spectral splits occur in the cooling 
phase. However, such splits occur at different energies for different post-bounce times 
and are smeared by multi-angle effects, disappearing in the DSNB by the integration 
over the time and over the SN population. 

• The dependence on the mass hierarchy is of the order 10 — 20% (Figs. 3 and 4) and 
it is strongest for antineutrinos. This is because, with the recently measured large 
value of #13, the high density MSW resonance is adiabatic and v e change from almost 
complete survival (survival probability ~ 0.7 — 0.8) to almost complete conversion 
(survival probability ~ 0.1 — 0.2) as the mass hierarchy changes from normal to inverted. 

• Combining results for different progenitor stars, instead of using the 10.8 Mq spectra 
for all stars, tends to increase the DSNB by ~ 5 — 10%, due to the more luminous fluxes 
of the more massive stars. 
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In summary, we find that, in first approximation, the DSNB can be described by a 
simplified model with a fixed, time independent spectrum for all supernovae, with MSW 
oscillations. The error due to neglecting the other effects that we have studied here is of 
the order of ten per cent or less and, assuming that the sign of the mass hierarchy will be 
known within the next decade and the picture of neutrino oscillations in supernovae will 
not drastically change in the next future, such error will be only related to the modeling of 
neutrino emission. We note that this error is smaller than astrophysical uncertainties: the 
error on the normalization of the supernova rate alone is at least 25%, and most likely higher 
once several systematic errors are included (among these, the still mysterious factor of ~ 2 
between the results of SN surveys and those of star formation measurements) . In addition to 
these, one should also consider errors at neutrino detectors, especially those due to the high 
level of background. 

In consideration of errors and experimental challenges, how likely is it that the effects 
we have discussed will have any relevance in future data analyses? When the DSNB will be 
detected, probably a first phase of data analysis will focus on excluding a number of models 
of neutrino spectra and SNR. In this respect, to model the DSNB as accurately as possible 
will be important to correctly establish or exclude compatibility. To actually test for effects 
at the 10% level will require long term progress to reduce at least some of the uncertainties 
of experimental and theoretical nature. For example, with new, precise measurements of the 
SNR and of the sign of the neutrino mass hierarchy from beams, and reliable neutrino spectra 
from future simulations (or from a galactic SN), it might be possible to use the DSNB to 
tests the luminosity of the neutrino fluxes from lesser known SN with high mass progenitors. 
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